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ABSTRACT 

Dust plays an important role in our understanding of the Universe, but it is not obvious yet how the dust in the distant 
universe was formed. I derived the dust yields per asymptotic giant branch (AGB) star and per supernova (SN) required 
to explain dust masses of galaxies at z = 6.3-7.5 (680-850 million years after the Big Bang) for which dust emission 
has been detected (HFLS3 at z = 6.34, ULAS J1120+0641 at z = 7.085, and A1689-zDl at z = 7.5), or unsuccessfully 
searched for. I found very high required yields, implying that AGB stars could not contribute substantially to dust 
production at these redshifts, and that SNe could explain these dust masses, but only if they do not destroy most of 
the dust they form (which is unlikely given the upper limits on the SN dust yields derived for galaxies where dust is 
not detected). This suggests that the grain growth in the interstellar medium is likely required at these early epochs. 

Key words, stars: AGB and post-AGB - supernovae: general - dust, extinction - galaxies: high-redshift - galaxies: ISM 
- quasars: general 


1. Introduction 

Dust plays an important role in our understanding of the 
Universe, both hiding the regions of intense star formation 
and allowing this information to be regained by the obser¬ 
vations of its thermal emission. It is, however, not obvious 
yet how the dust in the distant universe was formed, as 
dust formation requires specific conditions of relatively low 
temperature and high density. These conditions are met in 
atmospheres of asymptotic giant branch (AGB) stars and 
expell ed shells of supernova (SN) remnants (see I Gall et al.l 
for a review). Another possibility is that these stel¬ 
lar sources produce only dust seeds, and most dust mass 
accumulation happens in the interstellar medium (ISM; 
Draine fc Salpeteilll979D . 

Du st has been observe d to form in numerou s AGB 
stars (iMeixner et al. 200^ Matsuuraet^y_ 200^ 20^; 


a 


Sloane^_^ 2009|; ^inivasan et al.l[20(M iBover et al.112011 

20121 : iRiebel et ^ 2012fl . and it was shown theoreti¬ 

cally that one AGB star can p r oduce up to 0.04 

of _ du st (iMorgan fc Edmund: 

20061: IVentiira et al ' 


f r oduce 

200l 


20121 : Nanni et al 


Ferrarotti fc Gaill 

201.1 


120141: 


Schneider et al.l 2014 ). This maximum dust mass is re- 


Theoretical models predict t hat a SN can p r oduce 
at most ^ 1.3 Mq of dust (|Todini fc Ferrara 120011 : 

iNozawa et al.l [^00,1 . but likely only < 0.1 Mq survives 
the associated shock s and is released into th e ISM 
(iBianchi fc Schneider 20071: iGherchne ff fc Dwei3 120101 : 


iGall et aTT 2011al : iLakicevic et al.l I2015l1 .~ Large amounts 
of dust have b een fo u nd in the remnants of Cassiopeia A 

(|Morefan et al.l 

1987A (- 0.4- 


( Dunne et al.l 2?) 0^ 2009ll 


12001 iGomez et a1.li2009|) 


and Kepler 
IMa); SN 


0.7 mL: iMatsuura et al. I2011I : lindebetouw et al.l 20141) . 
and the Crab Nebula (~ 0.02-0.2 Mr^ iGomez et al. 120121 : 
iTemim fc DwekI 120131) . However, the high du st yields for 


Cassiopeia A and K e pler are still controv e rsial (iDwek 2004 • 


indKep _ _ __ 

Krause et al.l (200^_lGomez_et^r _ 200^_lWllsoi^_&_^atrlu 


^ Blair et al. boOTl: ISibthorpe et all 20101: iBarlow et ~ 

3) , and dust masses for other SN remnants are typically 

-3_in-2 


2004: Borkowski et al.l 20061: 

Sugerman et al. 

2006; 

Ercolano et al. 20071 Meikle et al. 

|2007 

: Rho et al. 

2008, 

20091: Kotak et a 

.1 20091: Lee et al. 

20091; 

Sakon et al.l 

2009; 

Sandstrom et al.l 

2009: iWesson et al. 120091: iGallagher et al. 

20121: Temim et al.l 20121). Some of this discrepancv can be 


turned only for a narrow range of the mass of an AGB 
star proge nitor (^ IMcA and for super-solar metallicity 
(Table 3 in iFerrarotti fc Gailll2006t) . For solar and subsolar 
metallicities, up to 0.02 Mq of dust can be produced by 
an AGB star, with a typical value of ~ 0.001 Mq across 
a wider mass range. It has been claimed that AGB stars 
produce most of the stellar dust (as opposed to dust grown 
in the ISM) in the Milky Way a nd local galaxies l|Gehr j 
ll989l : [^ ukovska^Henni^l2013). but also in high-redshift 
quasars ( Valiante et al.ll200a 20111) . 


Scottish Universities Physics Alliance 


attributed to different observed wavelengths used by differ¬ 
ent authors. If the dust is relatively cold, then mid-infrared 
Spitzer observations do not probe the dominant cold 
dust component, and would underestimate the total dust 
mass. Moreover, this difference may be connecte d with the 
evolut ionary stage of a SN remnant. Recently, iGall et ^ 
(|2014l1 found very large grains in the circumstellar medium 
of a SN just a few months after the explosion, and claimed 
that such grains are able to survive the shocks and help to 
condense larger amounts of dust a few tens of years later. 

It has been claimed that SNe are the main dust produc¬ 
ers at high redshifts, when at least some of the intermediate- 
mass stars had not had enough time to reach the AGB phase 
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Table 1. Properties of the galaxies in my sample. 


Name z log(Mdust/Af©) / Mq) log(Mdyn/M©) \og,(Mga_^ / Mq) Ref 


HFLS3 

ULASJl 120-1-0641 
A1689-ZD1 

6.340 

7.085 

7.500 

8 75+^'^^ 

‘ ^- 0.11 

■7 '7/-)-(-0.09 
'^- 0.11 

1 -OG.o.SO 

<7 q'7”I"0.09 

'-O'_ 0.11 

‘ •^^-0.48 

8 41+°-^'’ 

‘ -OO-O.IS 

10 57+^-^^ 
iU.O ( _o.30 

Q qo+0.05 

10.041-- 

9-3o1°:S9 

11 46 + 0-30a 
^-^•^'-'-0.30 

^•^^ — 0.09 

11.43 

11.00 

11 02'''^'^^ 
10.081“;-'’ 


HCM6A 

6.560 

< 7.44 



9 

^•*^^-0.30 





5,6 

Himiko 

6.595 

< 6.94 



iu.zy_Q Qg 

^■^^ — 0.02 




7^^ 

A1703-ZD1 

6.800 

< 7.19 








10.6 

IOK-1 

6.960 

< 7.27 



9 44+0-30 

^ —0.30 





11.6 

Z8-GND-5296 

7.508 

< 8.11 



9 44+0-30 

^ —0.30 





12 ^ 


Notes. Unless noted otherwise the values are from the references listed in the last column. St ellar masses are listed for the 
IChabrieil (120031 1 IMF. I calculated the stellar mass from the p hotometry (iRiechers et al.l[20131 1 using the Grasil code. I 
calculated the gas mass from the Lcii luminosity (IVenemans et al.1l2012ll using the empirical calibration o f ISwinbank et aO ll2012ll : 


Afgas = (10 ± 2l (Lcii/Lf7il. 


Reference s: 111: IRiech ers et a ll (l2013ll, H: Mortl ock et al ] (l201lll. B Venemans et al.l ll 2012 ll . 


Hu et al.l (l 200 ^. hi: ISchaerer et al 

Bradley et al.l (|2012ll . IllUlve et al 


7l: lOuchr e^^ I (l 200 ^. Isl: iHirashita et al.] (l2014ll . 
El: iFinkelstein et al! 


0: I Watson et al.l (l2015ll, H 
1: Ouchi et al.1 (1201, 3ll . flC 


(iDwek et al. 

2007: Dwek & Gherchneffl I 2 OIII: Dwe 

c et al .1 

20111 12014 

Michalowski et al.1 l2010cllhl: Gall et al.1 

2011 hi 

Hiorth et al. 

I20l4 iRowlands et al.1 1201411. However, us- 


ing the ALMA non-detection of a z ^ 6.6 Lya emitter 
(jOiichi et al.ll 2 ni^ . IHirashita et al.l (|2014r) found that only 
< 0.15-0.45 M 0 of dust can be returned to the ISM by a 
SN. 

Because of this difficulty in finding an efficient 
dust production mechanism by stellar sources, the 
grain growth in the ISM has been claimed to be nec¬ 
essary for the Mi l ky Way (iDraine &: Salpet^ 119791 : 
iDwek fc Scald 1 19801: ibrainel 199(1 ^09(1. loc^ galax¬ 
ies (|Matsuura et al.l 120091: Grootes et al.l 1201, ill, and 


distant ga laxies (iDwek et al.l 120071 : 


2010^2;__|Hirashit^&K^ |2011 


Michalowski et ahl 

MattesoJ 


Mattsson fc AndersenI l2012t jRuo&jlirashital 

Asano et al. 120131: ICalur a et al.1 [201 4 iHiorth et ^ 201A 

Rolands et al.| 2014 IValiante et al.1 I^l4 iNozawa et al.l 

201511 . The~ grain g rowth timesc ale _ is relatively short 

(a few tens of Myr; iDrainel ll990L 120091 : IHirashita | 2000 L 
I 2 OI 2 I : IZhukovska et al.1 l2008ll , which is consistent with 
very little variation in the metal-to-dust ratios in gamma- 
ray burst (GRB) and quasar absorbers, as this implies 
very rapid dust produ ction soon after metal enrichment 
(jZafar fc WatsonI 1201,ill . However, metal-to-dust ratios are 
higher at metalli cities well below 0.1 s olar, as shown by 
IZwl 8 (Fig. 2 in IZafar fc Wats^l201.11 and by a steeper 
slope of the gas-to-dust ratio for low metallicities (Fig. 4 in 
iRemv-Ruver et ^1201411 . This is connected with the grain 
growth timescale being much longer in th e metal-poor 
environments (Eq. (20) of lAsano et al.ll2013^ when stellar 
sources are the main dust producers. 

Apart from using the total dust budget arguments, these 
three possibilities are very difficult to distinguish because 
doing so would require precise measurements of dust prop¬ 
erties (grain size distribution and chemical composition). 
Therefore, SN-origin dust has only b een claimed based on 
extinction curves of a z ~ 6.2 quasar (iM aiolin o et al.ll^04 
iGallerani et al.1 12 OIOI : but see |Hiorthjet_^ 201^ and two 
GRB host galaxi es at z ^ 6.3 (IStra tta et aL[ 2007|: but see 
IZafar et al.il201ClLl2011h and z ^ 5 ( Perlev et al.ll201(ill . 


The objective of this paper is to investigate whether 
SNe and AGB stars are efficient enough to form dust at 
6.3 < z < 7.5 (680-850 million years after the Big Bang), or 


if grain growth in the ISM is required. I use a cosmological 
model with Hq = 70 km s“^ Mpc“^, Ha = 0.7, and = 
0.3. 


2. Sample 

For this study I selected all galaxies at z > 6.3 for 
which dust emission has been detected (with the excep¬ 
tion of SDSS J1148+5251 , a z = 6.42 quasar that was 
a subject of the study in iMichalowski et ^ l2010bL and 
MACS1149-JD, a z ~ 9.6 Lyman break galaxy [LBG], 
because its association with the 2 mm GI SMO object has 
not yet been confirmed: IDwek et al.ll20T^ . The sample in- 
cludes HFLS3, a re d Herschel-selecied z = 6.34 galaxy 
(|Riechers et al.ll201,‘^: IJLAS J1120+064T a colour-selected 
z = 7.085 quasar (jMortlock et al.l 120111 : IVenemans et al.1 
l2012f) : and Al689-z Dl, a lensed z = 7.5 Lyman break 
galaxy (jWatson et al.ll20T5h . 

Additionally, I selected galaxies at z > 6.3 for which the 
dust emission has been unsuc cessfully searched for: z = 6.56 
lensed Lya emitter HCM 6 A (|Hu et al.ll2002l : lKanekar et al.1 
I 20 TI magn ification factor of 4.5), z = 6.8 lensed LBG 
A1703-zDl (jBradlev et al.ll20T^ ISchaerer et al.ll2015L mag- 
nihcation factor o f 9), z = 6.595 Lya emitter Himiko 
(Ouchi^_etaL 20091 2 01^, z = 6.96 Lya emitter IOK-1 
( Ive et alT 20061: lOta et all 2014!) . a nd z = 7.508 LBG z 8 - 
GND-5296 (IFinkelstein et al.1 20]5 ISchaerer et al]l2015^ . 

Table [T] presents the relevant properties of these galax¬ 
ies, which are adopted from the references listed in the 
last column, except for the third stellar mass estimate for 
HFLS3 and the gas mass estimate for ULASJ1120+0641. 
Dust masses from ISchaerer et al.l (|2015!l are adopted with 
the Tdust = 35 K option. Multiple values for a given prop¬ 
erty denote different assumptions adopted when deriving 
it (temperatures and emissivity index for dust masses, and 
star formation histories for stellar masses). I converted all 
dust mass estimates to a common mass absorption coeffi¬ 
cient of Ki 2 oo/im = 0.67cm^g“^, and stellar masses to the 
Ghabrier IMF. 

For the stellar mass es timation for H F LS3 I used the 
photometry presented in IRiechers et al.1 ( 20131) a nd ap¬ 
plied the SEP fitting m ethod detailed in IMichalowski et all 
(l 2 (M I 20 M I 2 ninalld. see therein the discussion of the 
derivation of galaxy properties and typical uncertainties). 
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Fig. 1. Dust yield per AGB star (top) and per SN (bottom) required to explain dust masses measured in z > 6.3 galaxies for 
all possible combinations of dust and stellar masses (Table [T] the points were shifted slightly in redshift for clarity). Symbols are 
colour-coded according to the way the stellar mass was estimated: from spectral energy distribution modelling (black), assuming 
that the stellar mass is the difference between the dynamical and gas masses (yellow), or that it is equal to the dynamical mass 
(red). Filled circles, crosses, and open circles denote the assumption of the Chabrier, Salpeter (a = 2.35), and top-heavy (a = 1.5) 
initial mass functions, respectively. Small symbols with arrows correspond to galaxies with only upper limits for dust masses (lower 
part of Table [T]) . 


which is based on 35 000 te mplates from the library of 
Iglesias-Paramo et al.l ^^2007^ and some t e mplat es from 
Silva et al.l (|1998ll and iMichalowski et a hj (|2008[1. all of 
which were developed using GrasiiQ ( Silva et ^ I1998D . 


They are based on numerical calculations of radiative trans¬ 
fer within a galaxy, which is assumed to be a triaxial ax- 
isymmetric system with diffuse dust and dense molecular 
clouds, in which stars are born. 

I estimated the gas mass of ULASJ1120+06 41 from the 
[Cii] luminosity, Lcii (|Venenians_et_ah I2012D. using the 
empirical calibration of Swinbank et al.l 1 201211 : Mgas = 
(10 ± 2)(Lcii/.^0)- The uncertainty of this estimate has 
little effect on my result, as this mass corresponds only to 
^ 10% of the dynamical mass estimate. 

3. Method 

The redshifts of the galaxies considered here correspond to 
at most ~ 500 Myr of evolution, depending on the forma- 

^ http://adlibitum.oats.inaf.it/silva/grasil/grasil.html 


tion redshift. Hence, for AGB stars 1 only considered stars 
more massive than 3Mq, as less-massive stars have not 
had time to leave the main sequence and to produce dust 
(the main-sequence lifetime i s 10^° yr x e.g., 

iKipoenhahn fc WeigertlHo^ . 


1 calculated the dust yields per AGB star and per SN 
(amount of dust formed in ejecta of one star) required to 
explain th e inferred dust masses in these galaxies as de¬ 
scribed in IMichalowski et al.l (|2nincllbl ). Namely, the num¬ 
ber of stars with masses between Mq and Mi in the stellar 
population with a total stellar mass of Mgtar was calculated 
as N(Mo-Mi) = Q M-^dM/ M-^^MdM, 
where (Mq,Mi) = (3,8)M0 for AGB stars and (8,4O)M0 
for SNe, respectively. Depending on the galaxy, Mgtar was 
estimated in three ways: derived from the SED mod¬ 
elling (reported in Cols. 6-8 of Table [T]), assumed to be 
equal to Mdyn — A/gas and equal to Mdyn- The average 
dust yield per star is M(iust/A^(Afo~Afi). With these as- 
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sumptions 


this 


yield is 
MdM/ 


/ 


(Md ust /Mstar), where / = 
JMo (TableHI). 

I adopted an initial mass function (IMF) wit h Mmiu = 
0.15, Mma.y = 12 0 Mdi, and three possible shapes: fChabrieti 
l|2003ll . Salpeteil (1195^ with a slope a = 2.35, and top 
heavy with a slope a = 1.5. 


4. Results 

Figure [1] presents the dust yields per AGB star and per 
SN required to explain the observed dust masses using all 
possible combinations of stellar and dust masses presented 
in Tableland three choices of IMFs. 

The maximum theoretical yield for an AGB star is 
shown as 0.04 Mq because, as explained above, it is a strict 
upper limit. For SNe I plot the theoretical value of a dust 
yield with no dust destruction (1.3 Mq), and the expected 
dust mass returned to the ISM if dust destruction is taken 
into account (0.1 Mq). 

Among the three galaxies with dust detections Fig. [T] 
shows that only the dust mass of ULAS J1120+0641 can 
be explained by both AGB stars (required yields 0.015- 
0.04 Mq for the low dust mass options), or SNe with 
dust destruction taken into account (required yields 0.03- 
0.1 Mq), and only if the real dust yields per star are close to 
the maximum theoretically allowed values. However, AGB 
stars are not able to produce dust in two other galaxies 
(required yields 0.1-1 Mq), and SNe would need to return 
most of the formed dust without destroying it (required 
yields 0.2-1 Mq). 

Most of the galaxies not detected at the dust contin¬ 
uum (lower part of Table [T]) have stellar masses that are 
too small to provide meaningful constraints on AGB and 
SN dust yields, in the sense that the dust mass expected 
from the resulting number of these stars is much lower than 
the measured dust mass upper limits. In Fig.[Tl I show only 
the upper limits for Himiko and IOK-1, which are massive 
enough, so that the expected number of AGB stars and 
SNe is high enough to give reasonable constraints on the 
maximum dust yield per star. I found comparable results 
for Himiko with Hirashita et al.l ([203); in this galaxy each 
SN could have produced only < 0.2 Mq of dust assuming 
the top-heavy IMF, and < 0.03-0.09 Mq for the Chabrier 
and Salpeter IMFs. In IOK-1, one SN could have produced 
< 1.1 Mq of dust with the top-heavy IMF, and < 0.5 Mq 
with the Chabrier and Salpeter IMFs. The dust yield upper 
limits for AGB stars are not constraining (above the the¬ 
oretical limits), except for Himiko with the Chabrier and 
Salpeter IMFs, for which I derived < 0.01-0.03 Mq of dust 
per AGB star. 

Table [2] presents the factor / by which the Afdust/Mstar 
ratio needs to be multiplied for a given IMF to obtain the 
dust yield per star, and the factor by which the SED-derived 
stellar masses were multiplied for this IMF (the s tellar mass 
facto r for the top-heavy IMF with a = 1.5 is from lSternberd 
Il998ll . These factors can be quickly used to calculate the 
average dust yields, and they only depend on the adopted 
IMF and stellar mass ranges of dust producers. 


5. Discussion 

The dust yields per star shown in Fig. [T] imply that dust at 
z 6.3-7.5 could not have been entirely formed by AGB 


Table 2. Initial mass functions assumed. 


IMF 

/agb“ 

/sn“ 

Mstar/Mchab^ 

sym.'^ 

Chabrier 

29 

84 

1 

• 

Salpeter 

41 

127 

1.8 

X 

Top-heavy 

47 

54 

0.32 

o 


Notes. IMF-dependent factor, defined such that the re¬ 
quired dust yield per star is / x (Mdust/Mstar). The assumed 
stellar mass ranges are 3-8 Mq for AGB stars and 8-40 Mq for 
SNe. *■*’1 Ratio of the stellar mass derived using this IMF to the 
Chabrier IMF. Symbol used in Fig. [T] 


stars, and that SNe are efficient enough only if they do not 
destroy the pre-existing dust and most of the dust they 
produce. However, such high SN dust yields are difficult 
to reconcile with upper limits derived for galaxies not de- 
tect ed at the dust continu um (Himiko and IOK-1, see Fig.[T] 
and iHirashita et'aLll2014| l. This suggests that most of the 
dust mass in these galaxies had been formed in the ISM. 
This is in line with similar claims for galaxies at lower red- 
shifts. If the grain growth in the ISM is as fast as ~ 10 Myr, 
then even at these redshifts there is a sufficient amount of 
time for such dust mass accumulation (metallicities of the 
galaxies considered here are likely to be above the critical 
value of ~ 0.1-0.3 solar, the point at which the contribu¬ 
tion of the grain growth falls below that of st ellar sources 
because of the long timescale of this growth; lAsano et al.l 
[ml:lR.emv-Ruver et ahllMl . 

These conclusions obviously depend on how accurately 
the relevant masses were measured. I show here that they 
are robust because all my critical assumptions are conserva¬ 
tive in the sense that they tend to bias the resulting required 
dust yield per star towards low values: i) 1 assumed a rela¬ 
tively hig h value of the mass absorption coefficient k (com¬ 
pare with lAlton et al.l 1200411 . higher than what is usually 
assumed for distant galaxies, resulting in a systematically 
low Mdust; I included the stellar mass alternatives, which 
represent virtually the maximum value, and I assumed that 
the entire dynamical mass is composed of stars (red sym¬ 
bols in Fig.[T]) and assumed a double-component star forma¬ 
tion history (SFH) in the SED modelling, which has been 
shown to r esult in higher stellar masses th an single compo¬ 
nent SFHs (|Michalowski et ^l2012Ll2014^ : Hi) by including 
all stars more massive than 3 Mq (with the main-sequence 
lifetime of < 650 Myr), I implicitly assumed that the entire 
stellar mass was formed at the very beginning of the galaxy 
evolution, because otherwise a fraction of stars formed later 
had not finished their main-sequence stage and so should 
not be taken into account (the final dust mass depends on 
the SFHs, not only on the final stellar mass). This extreme 
SFH is not ruled out by the data, but is highly unlikely; iv) 
the theoretically allowed (coloured) regions in Fig. [1] rep¬ 
resent the maximum dust yield per star, which does not 
apply to all stars in a given mass range, so the average dust 
yield per star is lower; and v) 1 implicitly assumed that 
dust injected into the ISM is not subsequently destroyed 
by astration, interstellar shocks, or outflows, otherwise the 
produced dust masses would be higher than the currently 
existing dust masses and, in turn, the required stellar dust 
yields would be higher. In each of these cases different (more 
realistic) assumptions than those chosen here would result 
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in even higher required dust yields, making my conclusions 
even stronger. 

Moreover, in this analysis I considered AGB stars and 
SNe separately as dust producers. In reality, both sources 
operate simultaneously, at least when AGB star progenitors 
evolve off the main sequence. If the contribution of AGB 
stars and SNe is equal, then this would bring the points in 
Fig. [1] down by a factor of two, which would not change my 
conclusions. 

The choice of the IMF does not influence these conclu¬ 
sions strongly. For the same dust and stellar masses the 
Salpeter IMF results in yields that are a factor of ~ 1.5 
higher than for the Ghabrier IMF (see factor / in Table [5]), 
which makes it even more difficult to explain the observed 
dust masses. Since the SED-derived stellar masses are a 
factor of ^ 1.8 higher for the Salpeter IMF, the resulting 
yields using these SED-derived masses are almost identical 
to those derived with the Ghabrier IMF. 

Adopting the top-heavy IMF (and keeping masses the 
same) would increase the required dust yield per AGB star 
by a factor of 1.6 (/agb in Table [5]), virtually ruling out 
their contribution to the dust production in these galax¬ 
ies. On the other hand, in the top-heavy IMF the fraction 
of the most massive stars is higher, so the required dust 
yields per SN are lower by a factor of ^ 0.65 (/sn in Ta¬ 
ble [5]). However, the SED-derived stellar masses are 3 times 
lower, so all in all the number of these massive stars is lower 
than in the Ghabrier IMF, so the required dust yields are 
higher. Hence, the Ghabrier IMF is the most conservative 
assumption (resulting in the lowest required yields). 

Finally, it needs to be pointed out that the three galax¬ 
ies with detected dust emission are very different to each 
other, as the sample includes a dusty starburst, a quasar, 
and a more typical star-forming galaxy (the galaxies with 
dust non-detections also fall into this last category). This 
means that their dust properties and origins may be vastly 
different, even at similar stellar masses. Hence, the relative 
contribution of AGB stars, SNe, and the ISM grain growth 
may be different in these sources. 

6. Conclusions 

Using all three galaxies at z = 6.3-7.5 (680-850 million 
years after the Big Bang) for which dust emission has 
been detected (HFLS3 at z = 6.34, ULAS J1120+0641 at 
z = 7.085, and A1689-zDl at z = 7.5) and galaxies at the 
same redshift interval for which dust emission has not been 
detected, I found that AGB stars could not contribute sub¬ 
stantially to the dust production at these redshifts, and that 
SNe could explain their dust masses, but only if they do not 
destroy most of the dust they form . This last constraint is 
unlikely given the upper limits on the SN dust yields de¬ 
rived for galaxies where dust has not been detected. This 
suggests that grain growth in the interstellar medium is 
likely to be required at these early epochs of the evolution 
of the Universe. 
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